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Abstract. We are concerned with the question of determining the set C(Q), where C is a curve 
defined by an equation of the form y q — f(x), where q is an odd prime and / is a polynomial defined 
over Q. This question can often be answered using a set which encapsulates information about local 
solubility of a particular collection of covers of C. We define this set and show how to compute it. 



1. Introduction 

When working with the set of rational points C(Q) of an algebraic curve C, we often encounter 
examples where the Hasse principle fails to decide whether this set is empty or not. The method of 
two-cover descent on hyperelliptic curves is based on the fact that, for a hyperelliptic curve C, there 
is a computable collection of covers (j) a : Da — ^ such, that 

C(Q)= [J <pa(D a (Q)). (1.1) 

aGfinite set 

Therefore, when local-to-global arguments cannot be applied directly to C, the problem can be trans- 
ferred to the one of looking for rational points on the covers. This is described explicitly in [6j. In 
the present paper we extend this method to superelliptic curves C defined by an equation of the form 
U q = f( x ): where q is an odd prime and / £ Q[x] is q-th power-free. The theory behind the process 
of performing descent on the Jacobian variety Jq of such curves, defined over a field containing the 
relevant roots of unity, is studied in detail in [12J. In [7], an extension of the descent map to the Picard 
group of curves of this type is used to introduce new insights on the nature of LU( Jc/Q.)- Nevertheless, 
it is often much simpler and faster to avoid working with the rational points of the Jacobian and 
instead restrict to information obtained using only the initial curve. The Selmer set we define contains 
information which is sometimes sufficient to determine the set C(Q). We explain the necessary theory 
and present an explicit algorithm, similar to the one in [6], to compute this Selmer set. Note that we 
can not use a trivial extension of the existing routines because our algorithm is expected to deal with 
the possibility of singular points, which do not appear in the case of hyperelliptic curves. 

The problem of finding points on superelliptic curves generalizes the problem of finding solutions 
to Thue equations (see |10| ) and can also be used to solve generalized Fermat equations. In Example 
14.21 we solve four such equations considered in [9], which were used by the authors as examples of the 
limitations of their approach. Thus in many situations descent arguments are more appropriate than 
other techniques. As illustrated by Examples 14.11 and 14.21 it is often the case that C is everywhere 
locally soluble, but its associated covers fail to be so, preventing C from having any rational points, 
since the union in (|1.1|) is comprised of empty sets. In Example 14. 3| local information together with 
information obtained using subcovers, following a method proposed in [3] and [5], is used to prove that 
the curve 

y 3 = 0r 2 -3)(x 4 -2) 
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has no rational points except from one rational point at oo. By making descent applicable to singular 
superelliptic curves, we were able to prove in Theorem 14.51 that the only pair (a, 6) E Z> satisfying 

is (1,1). 

2. Descent on superelliptic curves 

2.1. Preliminaries. 

Definition 2.1. Let q be an odd prime and n be a positive integer. We define a superelliptic curve 

to be a plane curve C defined as the locus of the equation 

y q = f(x) = a n x n + . . . + aix + a , 
for some polynomial / with coefficients in Q. 

Remark 2.2. When q \ n we think of C as having a single point at infinity, otherwise C has q distinct 
points at infinity. 

Proposition 2.3. Every superelliptic curve is birational to a superelliptic curve satisfying an equation 
of the form y q = fix), where deg(f) = n and q \ n. 

Proof. Suppose C : {y') q = f'(x') is a superelliptic curve and that deg(/') = m = iq + j with i E Z>o 
and < j < q. Pick a £ <Q such that f'(a) ^ 0. Define the polynomial h by f'(x') = h(x' — a) and the 
polynomial / of degree n = q{i + 1) by f(x) = h(l/x)x q ^ l+1 \ Let C be the superelliptic curve defined 

by y q = f(x)- Then C is birational to C via the map (a/, y') i— > (^ x i^_ a , ^ x i- a y+i J ■ ^ 

In light of this proposition, we will assume that a n ^ and that q \ n throughout the rest of the paper. 
When this holds, we can also think of C as being the locus of the equation 

Y q = F(X, Z) = a n X n + ...+ aiXZ n - 1 + a Z n (2.1) 

in the weighted projective plane P 2 (l, -, 1), where the variables X and Z have weight 1 and Y has 
weight n/q. We will be using the affine and projective descriptions of C interchangeably. Without loss 
of generality, we may assume that / is a q-th power- free polynomial with integer coefficients. 

From now on let k denote one of Q, Q, Q p , Q p (where p can be any rational prime). Define J\ to 
be the set of monic, irreducible over k, polynomials such that 

f = a n h n \ 

where 1 < < q — 1 for all h E J 7 ^. Denote the degree of h by dh- Let be the semi-simple k- algebra 
k[t]/(g(t)), where 

9=U h 

Note that g is defined over Q and does not depend on k. We denote its degree by d. decomposes 
as a direct product of finite field extensions of k 

A k = H K h = H k[t]/(h). 

Denote by 9^ E the image of the generator t under the quotient map and by Ok the set 

@k = {0 h ■ he J" k }. 

For the following definition we will assume that a point (X,Y,Z) E C(k) is normalized such that if 
k = Q then X,Y,Z £Z with gcd(X, Z) = 1 and if k = % then X,Y,Z E Z p with either Z = 1 or 
X = 1, Z E pL v - We can always find such representations by scaling the points. 
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Definition 2.4. For k = Q and k = Q p define the component maps 5h '■ C(k) — > K^/K^ q 

_( (x-e h z)K* h q ax-e h z^o, 

6 h (X,Y,Z)-^ n ^F h ( XjZ) -i K *i otherwise) 
where is the two- variable polynomial with coefficients in Kh defined by (X — 9f l Z) nh Ff l (X, Z) 



F{X,Z) and n \J F h (X , Z)- 1 K* h q is defined to be the unique element v G K* h jK* q such that v nh = 
F h {X,Z)-\ 

Remark 2.5. (1) Note that v exists because the groups K^/K^ q have exponent q and gcd(g , ,n/ l ) = 1 

and it is unique since these groups are F g -vector spaces. 
(2) The component maps 5h are defined this way because for a point (X, Y, Z) G C(k) with F(X, Z) / 

0, we have that {X-6 h Z)K* h q = "^F h {X, Z)" 1 ^, since (X - 6 h Z) nh F h (X, Z) = F(X, Z) = Y q . 
Definition 2.6. Let 5 k : C(k) -> A*JA^ be 

In order to account for the fact that for any A G k* and (X, Y, Z) G C(k), (X, Y, Z) = (XX, \ n / q Y, \Z) G 
C(k), we quotient the codomain of 5\ by this action of scalars. So we define an action of k* on A^ by 

A W k = ( Aq * W k > 

where A G k* and (a^)^^ G This action descends to an action of k* /k* q on A^/A^ 9 . We denote 
by Al/k*A^ q the quotient of ^/yl^ 9 by this action. 

Definition 2.7. Define the descent map to be the composition 

C(k) A* k /At q Ai/k*A* k q 

where tt^ is just the projection to the quotient. 

2.2. The image of Sq. The image of 5q is contained in a finite subgroup A(q, S) of A^/A^. To see 
this let us restrict our attention to finding the allowed possibilities for each of the jfcJ-Q components. 
Let h G J-q and set f h (x) = f{x)/(x - 9 h ) nh G K h [x]. 

Suppose (X, Y, Z) G C(Q) with X, Y, Z G Z and X coprime with Z. We have that 

4(x,y,z) = {x-e h z)K* h q . 

Now suppose that p \ a n OK h is a prime ideal of the ring of integers Ox h of the number field Kh- By 
assumption we have that 

ord p (X - 6 h Z) > 0, oi& p {F h {X,Z)) > 0. 

At this point we want to figure out which primes p appear in the factorization of (X — 6hZ), but not as 
a q-th. power. So suppose q \ ord p (X — O^Z). Since gcd(n^, q) = 1, this implies that q \ ord p (F^(X, Z)). 
In particular we have X = Q^Z and Fh(X, Z) = modulo p, which together give that Fh(0hZ, Z) = 
%n-n h f h (0 h ~j = o modulo p. But if Z = modulo p then also X = modulo p which contradicts 
coprimality, so we have that p G Supp(fh{0h)OK h )- By dropping the initial condition on p we have 
that if q \ ord p (X - 6 h Z) then p G Supp(a n f h (9 h )OK h ) Q Supp(A£>^J, where A = a n Disc(g). In 
other words 

(x-e h z)o Kh =p e 1 1 ...p?x q 

where {pi, . . . , p^} = Supp (a n fh(0h)OK h ^ , (ei, . . . , e{) G and X is a fractional ideal of Kh- 

Now define the sets of primes Sh = Supp(a n fh(0h)OK h ) for h G Fq. Then by the discussion above 
Image(4) C K h (q,S h ) where 

Kh(q, S h ) ■= {atK* h q : q | ord p (a) for all p ^ S h }. 
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Note that Kh(q, Sh) is a finite subgroup of K^/K^ q (a proof of this can be found within the proof of 
Proposition VIII 1.6. in [15J). Therefore we have that 

Imaged) C JJ K h (q, S h ) =: A(q, S) (2.2) 

which is a finite subgroup of Aq/Aq. 
2.3. The image of S^. 

Definition 2.8. Define the weighted norm homomorphism N^/ k : A k — > k* as 



k 



Since N^/^A^ 1 ) is a subgroup of k* q we also get a homomorphism iV^/^ : A^/A^ — > k*/k* q . 
By commutativity of the following diagram of norm homomorphisms 

.„ N A/k . . 

At k* 



At/A? ^ k*/k* q 



and the fact that 

Y q 

N A/k (S k (X, Y, Z)) = H N Kh/k (X - h Z) n » =Z n x]J h(X/Z) n » 



we can deduce that 



Image(4) C N^ k ( — k j : 11 



q 



which, if non-empty, is a coset of the subgroup Kerne^iV^/jj) in A^/A^. Combining this with the 
inclusion (12.21) we deduce that 



QH Q nA{q,S) =:H Q (S). 
2.4. The image of (iq. 

Lemma 2.9. Let Aq be the semi-simple Q-algebra associated to the curve C and A(q, S) be the subgroup 
of Aq/A*^ defined in ([2~3]) . Also set 

T := {p prime : for all h £ Tq, and all p £ Supp(pOK h ), (q | ord p (p0^ h )) or (p £ Sh)}- 

^ A* Q /A*^ 3 A* Q /Q*A% — > 1 
t A 

T) A(q,S) 



Then 

(i) Image(t) n A(q, S) = i (Q(q, T)) 
(it) n Q (A(q,S))^A(q,S)/t(Q(q,T)). 

Proof. (i) First suppose that aA* q £ i(Q(q,T)). Then q | ord p (a) for all p £T. Let p be a prime of 
Ox h for some h, with p ^ Sh- We know that 

oid p (aO Kh ) = oid p (pO Kh ) x ord p (a) 

so by the definition of T at least one of the factors will be divisible by q thus also their product, 
which implies that aA* q £ A(q, S). 

For the opposite inclusion, suppose that a £ Q* and aA* q £ A(q, S). Then q | ord p (aO for all 



Descent on superelliptic curves 



5 



h G and for all p ^ 5^ so if p ^ T, then g must divide ord p (a) since q divides the product but 
not the first factor in the equality above, 
(ii) We have 



Kernel(vrQ) n S) Image (t) D ^4(<?, S) t(Q(qr,T)) 



□ 



Let us denote ttq(Hq(S)) by JEZq(S). Since Image (<5q) C Hq(S), we have that 

Imaged) C Hq(S). (2.3) 

In the following section we will see how Image(/z<Q)) is contained in a potentially strict subset of Hq(S) 
and provide an algorithm to compute it. 

3. The Selmer set 

3.1. Determining the image of [iq . In this section we will provide an algorithm which determines 
Image(/XQ p ) for a rational prime p. The algorithm relies on the fact that points of C which lie in a 
"sufficiently small" p-adic neighborhood, have the same image under /Uq . 

The diagram below is crucial in the process of refining the possible image of \x even further: 

A* n /Q*A* q 



MQp > a* /o* A* q 

By commutativity, if we have a rational point P G C(Q) then o i p (P) = r p o ^q(P). Therefore we 
have that 

Image (//q) C r" 1 (Image (/xq p )) n Hq(S) (3.1) 

Definition 3.1. The Selmer set over Q of the superelliptic curve C, is defined as 

Sel^(C, Q) = {[a] G Hq(S) : r p ([a]) G Image(/ZQ p ) for all rational primes . 

Remark 3.2. Strictly speaking, the set defined here corresponds to the "fake " Selmer set found generally 
in the literature (e.g. [6] and [7]). Roughly, the difference between the fake and the actual Selmer set 
is that the latter distinguishes between covers (defined in Section T3.2p <p a : D a — > C and <f> a i : D a i — > C 
when (j) a and (j) a i are different even if D a and D a i are isomorphic. Since we are not using both sets, 
we omit the "fake" from the notation. 

After considering the inclusion (|2.3p in the end of Section 12.41 and the inclusions (|3.ip for every rational 
prime p we get that 

Image(/x Q ) CSel^(C,Q). 

Let h G J-q and denote by p^ the prime of Ox h - The following two lemmas are used to show that 
the analytic space C(Q P ) can be covered by a finite number of neighborhoods, where the map /j,q p is 
constant. In practice (or X) will be a finite precision approximation to the first coordinate of a 
point (X>,Y',l) G C(Q P ). 

Lemma 3.3. Suppose that X', Xj, G Z p with oid p (X' — Xf.) > k 

(i) ifk> 2ord p fa (g)+^ d pJ^-^)+ 1 then ( Xk - 9 h )K* h q = {X' - 9 h )K* h q . 

(ii) Ifk> 2°rdp h ( g )+°rdp,(7 h (^))+l ^ J h{Xk)K *« = f h ( X ')K* h q . 
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Proof. (i) By the assumption X' = X k + up k where u G Z p . So 

x'-e h np fe 



Now let r(t) = ~ t q . By Hensel's lemma we have that the following is a sufficient condition 



for r to have a solution in Kh, 



ord p ,(r(l))>2ord p , 



koTd Ph (p) - ord Ph (X k - 9 h ) > 2ord Ph (g) + 1 
So as long as the condition of the lemma is satisfied Hensel's lemma ensures that [X k — Oh) and 
(X' — 9h) are the same modulo K^ q . 
(ii) This is very similar to the previous part. Just use the fact that 

„fc\ _ A I Y. \ _L o,^ nrViara »i C ^ »U qo+ — fh( X ') 



fh(X') = f h {X k + up k ) = f h (X k ) + vp k where v £ Z p and set r(t) = - i 9 . 



□ 



Lemma 3.4. // {X k } k *L l C Z p is a sequence satisfying ord p (X' — X k ) > k for some X' G Z p and every 
k, then there exists N G Z>o such that Xn satisfies at least one of conditions (i) or (ii) of Lemma \3.3[ 

Proof. Suppose such iV does not exist. This means that for every k we have 

• ( 2ord Pfa (g) + ord p ,(X fc -^) + l 2ord Pfe (g) + oid Ph (f h (X k )) + 1 \ 

{ e ph/p e ph/p J 

and therefore both oid Vh (X k — 8^) and ord Ph (fh(X k )) tend to infinity as k tends to infinity. But since 
{Xk} r k L\ converges to X' we have that {X 1 — 6^) = fh(X') = 0, a contradiction. □ 

A clear distinction between the case of hyperelliptic (q = 2) and superelliptic (q > 2) curves is 
that, a superelliptic curve is allowed to have singularities, since / being q-th. power-free is no longer 
equivalent to / not having repeated roots. At this point we would like to use some version of Hensel's 
Lemma to determine whether our finite precision X k lifts to Z p as the first coordinate of a point 
(X' ,Y' , 1) G C(Q p ). We have to be careful not to ask this question for points approximating one of 
the singularities as that would result in an infinite loop. Thus we have to determine the size of the 
Hq p -constant neighborhood around each singularity (which is defined over Q p ) in advance and compute 
its image. 

Let 

= {h G Fq p ■■ deg(/i) = l,n h > 1, 9 h G Z p } , 
where the exponent Isi stands for linear, singular and integral. Elements of this set correspond to the 
singular points on C that are defined over Q p , but are of the form (9h, 0, 1) with 9h G Z p . The last 
condition arises because we split the computation into two parts, the first being the determination of 
the image under fj,Q p of points of the form (X 1 , Y', 1) G C(Q P ), with X',Y' G Z p . 
Consider the following functions: 

l: function SizeOfNeighborhood(/i) 

2: k h <- 

3: Finish <— false 

4: while [ Finish = false] do 

5: kh^kh + l 

6: X <r- X G Z C Z p : ord Ph (X - 8 h ) > k h 

7: VALLIST^- { 20rd ^' (g)+ ;; fe d ;;; (X -^ )+1 : h> G F Qp \ {h}}u{2ord Ph ( q ) + ord Ph (f h (X)) + l} 

8: if [ max(VALLlST) < k h } then 

9: Finish <— true 
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10: end if 

11: end while 
12: return kh, X 
13: end function 

Note that the function SiZEOFNEIGHBORHOOD only makes sense when deg(h) = 1 otherwise we 
would not be able to find an X satisfying the condition of step 6 for every given kh (that would imply 
that Oh £ and we will actually only apply it to elements of J-q 1 - This function has a double use: 
The returned value of X will be used to compute fJ,Q p (0h, 0, 1), and kh will keep track of the size of the 
fiQ p -constant neighborhood around the singularity. Now for h G Tq* set 

U h = {X €Z P : ord p (X - 6 h ) > k h } 

and 

The following function can be thought of as partitioning Z p into neighborhoods, with the partition 
becoming finer close to the singularities. Then the function LoCALlMAGE will test each of these 
neighborhoods for elements that lift to points on C, and if necessary partition them further into 
fj,Q -constant parts. 

l: function ComputeInputList(Reps, k, List) 
2: for X e Reps do 

3: if [ 3h € Jgj :[k h >k and ovd p (X - h ) > k h ] ] then 

4: List ^ComputeInputList({X + tp k : t e {0, . . . ,p - 1}} , k + 1, List) 

5: else 

6: List ^— List u{(X,k)} 

7: end if 

8: end for 
9: return List 
10: end function 

At this point we should stress that the functions SizeOfNeighborhood and ComputeInputList 
would be redundant if there were no singularities, and the function LoCALlMAGE would be sufficient 
to compute the local image. In that case the input (LlST={0, . . . ,p — 1},IMAGE={}) would produce 
the require result. 

l: function LocalImage(List, Image) 

2: for [ (X, k) £ List ] do 

3: if [ X g U or $(X', Y', 1) G C(Q P ) : ordp(X' - X) > k ] then 

4: List <- List\{(X, k)\ 

5: else 

' ' ' 2ord Vh (q)+ord Ph (X-e h )+l 2 ord Pfe (g)+ord Pfe (f h (X))+l 
e Ph/P e t>h/P 

7: if [ V[fci,A;2] GValList mm(ki,k 2 ) < k ] then 

NewElement ^— 1 
for [ [ki,k 2 ] GValList ] do 
if [ k\ < k 2 ] then 

NewElement <— NewElement x (X - 9 h )K* h q 
else 

NewElement^— NewElement x n \J ' fh(X)- 1 K* h q 
end if 
end for 

NewElement^— ttq p (NewElement) 



6: VALLlST^— 



9 
10 
11 
12 

13 
14 
15 
16 



he Tin 
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17: Image <— Image u{NewElement} 

18: else 

19: NewList^ {(X + tp k ,k + l) :t€ {0,...,p-l}} 

20: Imaged— LocalImage(NewList, Image) 

21: end if 

22: end if 

23: end for 

24: return Image 

25: end function 

With the help of the SizeOfNeighborhood function, we pre-compute and store in the variable 
image the images of the singular points under ^Q p . Using LoCALlMAGE with initial input 

(ComputeInputList ({0, . .. ,p— 1}, 1, {}) , image) 

we get as output a set V\ C Aq /Q*Aq which satisfies fJ,Q p (U{) = V±, where U\ = {(X, Y, 1) G C(Q P ) : 
1,7 6 Z p }. With slight modifications to the routine above we can also obtain as output a set Vi such 
that fJ.Q p (U2) = V2, where U2 = {(l,Y,pZ) G C(Q P ) : Y, Z £ Z p }. Thus we obtain the complete image 
of fiQ p , since U\ U U 2 = C(Q P ). 

Remark 3.5. We can be certain that the routine Locallmage terminates after a finite number of steps 
because of Lemma [3.41 An infinite loop would correspond to a sequence {^fclfc^i converging to some 
X' G Z p satisfying X' — 8^ = fh(X') = which is impossible. Also note that the If statement at step 
3, ensures that if X G U, in other words if our approximate value is very close to the first coordinate 
of one of the singular points, then it is excluded from LIST and we do not try to lift it using Hensel's 
Lemma, which would have resulted in an infinite loop. 

3.2. The corresponding covers. For every a G A^ such that [a] := aQ*A* q G we can construct 
an unramified cover of C of degree q d ~ 2 



:D a ^C, 



defined over k satisfying the properties 



D a (k)^0 [a] € Image (/%) 

[a] = [a 1 ] => D a = D a r 

First let us give an equivalent description of the k-algebra A k - Denote the absolute Galois group 
Gal(k/k) by Q k . We fix embeddings k for all h G J\ that are compatible in the sense that they 

agree on the intersections n K^i for h,h' G J\. We then get an inclusion G^ «— )■ Gg and we can 
treat elements at G as elements of k and elements of Qk as elements of 0^. 

Lemma 3.6. 

A k ^Ma Pk (0 E ,k), 

where the right hand side is the set of all Qk-equivariant maps from 0j. to k. 
Proof. The isomorphism is given by 

{ah)heT k ^ (Oh' ^ a OL h :a£G k , h G T k , a 6 h = 6 h >) 

with inverse 

^(^(0h))e h ee k ce t - 

□ 
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Let a G such that [a] G H^. Since aA* q G H^, there exists v G k* with an^YA/kCaO = v 9 . Let -D Q 
be the variety in P d_1 x C defined by 

(KW 6 , (X, Y, Z)) G £ Q «*3A ^ s.t. Aa(0 h )u£ = X - ft Z (3.2) 

for all h & J-g and 

A n / 9 v JJ u\ h = Y 



We use a description of the covers D a similar, at least in terms of their ambient space, to the one found 
in [7] and |14) . We equip the first factor, P d_1 , with the twisted ^-action which permutes coordinates 
in the same way it permutes Og. In other words if a G C% satisfies a 6h = Oh', then a Uh = Uh>- It is 
then obvious from the definition that D a is actually defined over k. Projection to the second factor 
gives rise to the required covering map 

(j) a : D a -t C. 

Lemma 3.7. The map (f> a '. D a — y C is unramified of degree q d 2 . In particular D a is a curve. 



Proof. Let (X, Y, Z) G C(k) and suppose X — 8h Z ^ for some ho G J-g. Then Uh ^ so we can set 



X—Q Z 

Uh = 1- By doing this we also fixed A = — (d ■ So for all Uh with h ^ ho we have 



q = a(6 h0 )(X - 9 h Z) 

Uh a (9 h )(x - e ho zy 

Since we are over k, if there does not exist h G J\ such that X — Q^Z = 0, then there are exactly q 
different choices for the value of each of the d — 1 u^s. Once d — 2 of them have been chosen, the 
remaining one is decided by the relation 



The fact that there is a unique choice for the value of the remaining coordinate uses that gcd(n/j, q) = 1 
for every /i£ J|. On the other hand if X — OyZ = for some h' G J-g \ {ho} then = and there 
are q choices for the remaining d — 2 u^s. The extra relation in this case does not decide the value for 
any of them. We see that in both cases the fiber of <f> a over P G C(k) contains exactly q d ~ 2 points. □ 



Proposition 3.8. D a (k) ^ if and only if [a] G Image^k). Furthermore, if [a] = [a'] then D a = D a 

1*9 
1 k • 



over k. In other words, up to ^-isomorphism, D a only depends on the class [a] in A^/k*A^ q 



Proof. Let P = [(u h ) he ^, (X, Y, Z)) G D a (k). Projecting gives <j> a {P) = (X,Y,Z) G C(k) and if 
X - 9 h Z ^ V h G then 

4 (MP)) = ((X - h Z)K^) 6he = (\a h ulK* h «) dhe 
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therefore /%(0 a (P)) = [a]. On the other hand if there exists h' £ such that X — Oh'Z = then by 
definition of 5h> we have 




' l h' K* q 



Xa h tKp. 



So again /Lik(</> a (P)) = [a], since all the other components of 5k are evaluated without using cofactors. 

For the other implication, suppose [a] £ Image(/%). This means there exist X,Y,Z G Ok, A € I* 
and (3 £ Al such that 

\ a (e h )p(e h y = x-e h z 

for 8h £ Bfc. After conjugating these relations by elements of C% and using the fact that a and j3 are C%- 
equivariant, we obtain the corresponding relations for 9h £ ©£ \ 6k- Then (^(f3(6h))heJ r ^, {X, Y, Z)) G 
D a (k). Note that if for some /i' G Pk. we have X — O^Z = then d/,/ = 1 and the corresponding 
coordinate is equal to zero. 

For the last statement suppose that we have a, a' G At with [a] = [a'}. This implies that there 
exist A G k* and (5 G At such that a = Xa'/3 q . By definition of the covers it is not hard to see that 
we can then map D a to D a > via ((u h ) heJ r-, (X, Y, Z)) H> ((/3(9 h )u h ) h ^, (X, Y, Z)) which is clearly an 
isomorphism and it is defined over k since it is invariant under the twisted ^-action. □ 

Corollary 3.9. Let % be any subset of Hq(S) containing Image(^Qi) ; then 

C(Q)= (J 4>a (A*(Q)) . 

[a]eH 

In particular the above equality holds for % = Sel^ (C, Q) . 

Proof. Since for all a G Aq such that [a] G Hq, D a and <p a are defined over Q, we know that the 
right hand side is contained in C(Q). Also, from the proof of Proposition 13.81 we deduce that if we 
have P G C(Q) and //q(P) = [a] then there exists Q G D Q (Q) with <p a {Q) = P- On the other hand, 
if [a] G H \ Image (/xq) then D a (Q) = 0. □ 

Proposition 3.10. T/ie curves D a are non-singular. 

Proof. To show this let us restrict to the affine patch where Uh ^ and Z ^ for some ho 6 J|. We 
thus assume that Uh = Z = 1, label the elements of J-g \ {ho} using an index i £ {1, . . . ,d — 1} and 
rename m := u^, 0{ := and m := n^. to simplify the notation. The defining equations for D a in 
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this patch become 

a(pi){X - 0„K =a(9 )(X - 0$ for i e {1, . . . , d - 1} 

d-l 

V (x - e ) n / q n =y. 

1=1 

We get the following (d + 1) x (f matrix of partial derivatives which represents a linear map whose 
cokernel is the cotangent space of D a at a generic point of the affine patch. 

/ f \ t— r \ 

_ //i \ _ £) .//IS _^/i\_a _//i\ \..0 _ //i \ .. / n \ / /i \n/o-l I I n, 1 



a{6 x )u\ - a(9 ) a{9 2 )u q 2 ~ a(9 ) . . . - a(0„) » (f ) (* - 9 ) n/q - 1 J] <' 

i=l 

ga(fli)(jr-fl K -1 ■•■ uniOY-flo)" 7 '*? 1-1 II 



o ^(ftOpr-flojuT 1 ... o unaCx-^o)" 7 ^ 2 " 1 n <' 

i/0,2 

... 9 a(^_i)(X-ff )«5li tm d _i(X-0o) n/9 <^ 1_1 JJ 

i^0,d-l 

V ... -1 / 

(3.3) 

This matrix has rank d at every point of the affine patch. To see this note that the first d — 1 entries 
of the first row, can never be zero, since this would contradict the fact that the roots of g are distinct. 
For the same reason, at most one row can be identically zero and when this happens the d x d matrix 
obtained by deleting that row has rank d. A similar argument for all of the 2d affine patches covering 
D a shows that D a is non-singular. □ 

Proposition 3.11. If p is a rational prime such that p \ qA and a G Aq such that [a] S Hq(S), then 
Dp has good reduction, where f3 € Aq^ is such that \fl] = r p ([a]). 

Proof Since p does not divide q or the discriminant of g, a S Aq can be chosen such that ord p (a(^/ l )) = 
for every h £ and every p | p. Take j3 to be the image of this a under the inclusion (8>Q P : Aq —> 
Aq . Then [f3] = r p ([a]) and also ordp (f3(0h)) = for every h € and every p | p. Thus the defining 
equation 

\p{e h )u q h = x- e h z 

of Dp can be reduced to an equation modulo ph for every h £ J~q , where ph is the prime of Kh above 

p. We thus get an unramified cover <f>n : Dp — > C defined over ¥ p . Furthermore we know that Dp is 
non-singular since the reduction of the matrix of partial derivatives in fj3-3|) has full rank. □ 

Proposition 3.12. The genus G of the covers D a is equal to q d ~ 2 ~ <?) + L 

Proof. Let tp : C — > P 1 be the map (X,Y,Z) i-> (X,Z). This has degree q and is ramified above 
the d points (X,Z) satisfying F(X,Z) = 0. Since <f> a is unramified of degree q d ~ 2 , the composition 
1/3 o (f> a : D a — > P 1 has degree q d_1 and is ramified at the same points as ip. Each ramification point 
has exactly q d ~ 2 preimages, each of them with ramification index q. The Riemann-Hurwitz formula 
applied to ip o (p a yields 

2G , -2=(2Genus(P 1 )-2)deg(>o0 a )+ ^ (eg - 1) 

PeP 1 ,Qe(V-o0 cl )-i(p) 

2G-2 = -2q d ~ 1 +dq d - 2 (q- 1) 
G=q^{^^-q) + l 

□ 
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Definition 3.13. Define the set of useful primes to be 

S U n := \p rational prime : p I qA or -Jp-\ < 2G 

{ Vp 

Proposition 3.14. Suppose that p is a rational prime and p ^ S up . Then Hq(S) C r" 1 (image (//<q p )) . 

Proof. Let a G Aq such that [a] G Hq(S) and /3 G Aq^ such that = r p ([a]). By Propositions 
13.111 and 13,121 we know that Dp is a non-singular curve of genus G over ¥ p . Then by the Hasse-Weil 
inequality we have that #Dp(JF p ) > and we can use Hensel's Lemma to lift to a point in Dp{Q^ p ). 
By Proposition 13.81 this is equivalent to [f3] = r p ([a]) G Image(^Q p ). □ 

Corollary 3.15. 

Sel M (C,Q) = {a G Hq : r p (a) G linage^) /or a// p G S up } 

3.3. Computational efficiency. The groups ^(g, S), Q(q,T) and the homomorphism l defined in 
Sections 12.21 and 12.41 can be computed using commands implemented by Claus Fieker in the MAGMA 
computer algebra system |3j, so we can compute ttq (Hq(S)) = Hq(S). The bottleneck of the compu- 
tation is computing the class and unit groups of the number fields K^, for h G J~q, which are needed 
for the construction of A(q, S). 

Also although Corollary 13.151 indicates that the algorithm computes Sel^ (C, Q) in a finite amount 
of time, the size of S up is prohibitively large and in general we can only hope to get information using 
small primes. Nevertheless we can still put the algorithm in good use as in most cases bigger primes 
do not have a contribution in cutting down the set we already have. 



4. Examples 

In this section we give examples of how descent on superelliptic curves can be used to tackle some 
interesting number theoretical problems. Example 14. H is a preparatory example to demonstrate how the 
results from the local computations are obtained and combined to prove statements regarding the sets of 
rational points of the curve, or curves in question. In Example 14. 21 we show how descent can sometimes 
be the appropriate technique for solving generalized Fermat equations. In Example 14.31 we consider a 
superelliptic curve, which although covers a plane cubic curve in an obvious way {X \— > X 2 , Z i— > Z 2 ), 
the set of rational points of this curve is infinite, so it is impossible to construct the set of rational 
points of the superelliptic curve by pulling back rational points of the cubic curve. After using the 
algorithm described above we exclude all but one of the covers D a due to local insolubility and we 
manage to compute all the rational points of this remaining cover D\ by other means (found in [1],[5])- 
Finally, we compute C(Q) since C(Q) = 0i(Di(Q)). The reasoning of this example, i.e. using covering 
techniques to transfer the problem to a different type of curves, is also used in the proof of Theorem 14.51 
but unlike Example 14.31 an the curves involved can be defined over Q. Theorem 14.51 also demonstrates 
the significance of extending descent arguments to singular superelliptic curves. 

Example 4.1. Consider the curve defined by the equation 

y 5 = 2x 5 + x 4 + 2x 3 + x 2 + 3a; + 3. 

This is ELS, but after applying the algorithm to this curve we obtain the results shown in Table [TJ 



p 


1 


2 


3 


5 ... 17 


19 . . . 37 


41 


#H Q (S) fl rf^CCCQ,))) 

l prime 
l<p 


25 


25 


25 


2 ... 2 


1 ... 1 






Table 1. Descent computation for C : y 5 = 2x 5 + x 4 + 2x 3 + x 2 + 3x + 3. 
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Therefore 

H Q (S) n r 5 - x Guq 6 (C(Q 5 ))) n r^ 1 H 19 (C(Q 19 ))) D r^ 1 (^ 41 (C(Q 4 i))) = 
which proves that C(Q) = 0. 

Example 4.2. In jS], Halberstadt and Kraus, consider the following four generalized Fermat equations 

16a 7 + 87b 7 + 625c 7 = (4.1) 

11a 5 + 296 5 + 81c 5 = (4.2) 

27a 5 + 166 5 + 2209c 5 = (4.3) 

32a 7 + 816 7 + 187c 7 = (4.4) 

These have solutions everywhere locally, but appear to have no rational points. The authors explain 
how the modular approach fails to show that the set of rational points is empty. We show how one 
can use descent to tackle all four of them. 

Observe that the problem can be easily transferred to the one of finding rational points on superel- 
liptic curves. 



(1441) O (-6, 2a, -c) G Ci(Q), where G\ : Y 7 = 8(87 X 7 + 625Z 7 ) 
M & (-a, 3c, -6) G C 2 (Q), where C 2 : Y 5 = 3(11X 5 + 29Z 5 ) 
g3D & (~a, 2b, -c) G C 3 (Q), where C7 3 : Y 5 = 2(27X 5 + 2209Z 5 ) 
g3D (-&, 2a, -c) G C 4 (Q), where C 4 :Y 7 = 4(81X 7 + 187Z 7 ) 
Table [2] contains the results obtained when we perform descent on these curves. 



p 


1 


2 


3 


5 


7 ... 23 


29 


i prime 
Z<p 


49 













#F Q (S) p| rf^CCfeCQ,))) 

Z prime 
Z<p 















#F Q (S) p rfH^C^CQi))) 

Z prime 
Z<p 


5 


5 


5 


1 


1 ... 1 





## Q (S) p r / - 1 (/i Qi (C 4 (Qi))) 

Z prime 
l<p 


7 














TABLE 2. Results for the generalized Fermat curves 



Therefore d(Q) = for % = 1, 2, 3, 4. 

Example 4.3. Consider the curve C in P 2 (l,2, 1), defined by 

(X, Y,Z) £C ^Y 3 = (X 2 - 3Z 2 )(X 4 - 2Z 4 ). 

This curve has three points at infinity (1, 1,0), (l,p, 0) and (l,p 2 ,0), where p is a primitive cube root 
of unity. So we know that it has at least one rational point. After applying the algorithm to this curve 
we obtain the results shown in Table [3j 

The element of Hq(S) remaining is the image of the point (1, 1,0) under /mq which is equal to the 
identity element 1Q*^4* 3 . This corresponds to a cover fa : D\ — > C, where D\ is a curve in P 5 x C 
defined as in f|3 . 2[> . whose set of rational points is non-empty. We fix embeddings of Kr x 2_$\ and 
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p 


1 


2 


3 


5 


7 


11 


13 


17 ... 


#s Q (s) p| rrVozCcm))) 


243 


243 


9 


3 


3 


3 


3 


1 ... 


i prime 



















Table 3. Results for C : y A = (x 2 - 3)(x 4 - 2). 



and 



6 



Q and index the six elements of 0q as 0i = \/3, i?2 = ~~ \/3> $3 
-iy2. D\ is denned by the following relations 

((tti, ...,n 6 ), € £>i 44>3A / such that Xu 3 - = 

for 1 < j < 6 and 



^,0 4 = -^2,0 5 =i^ 



0,Z 



3=1 

This covers a curve E' of genus 1 in P 2 , defined over the number field L : = Q(03), given by the equation 

(U,V,W)eE'& V 3 = (U -$zW){U 2 + $lW 2 ). 
We have the following commutative diagram 

D X C ((u 1 ,...,u 6 ),(X,Y,Z)) ] >(X,Y,Z) 



E' 



(Xui (X - $ 3 Z)u 5 u 6 , Zui) l — > (X, Z) 



Note that k is a well defined rational map between two non-singular curves so it is actually a morphism. 
We can put E' into Weierstrass form via a linear transformation, by moving the point (03,0, 1) to oo. 
We obtain the Weierstrass model 

E : y 2 z — Mlyz 2 = x 3 - 64z 3 . 
The isomorphism of the two models is given by 

i> : E -> E', i>(x, y, z) = {ti 3 y, 20 3 x, y - M 2 z). 
Using the package MAGMA we find that the Mordell-Weil rank of E(L) is 1. Since 

we are not interested in all the L-rational points of E', only (U, V, W) E E'(L) such that U/W G Q. 
Determining these will give us C(Q). To solve this problem we can use "Elliptic curve Chabauty" [3]. 
Fortunately this is implemented in MAGMA. Using the inbuilt MAGMA commands we find that 

{(U,V,W) G E'(L) : U/W GQ} = {(1, 1, 0), (0, — 3 , 1)}. 

We deduce that C(Q) = {(1, 1,0)}. 

Example 4.4. In [2], [5], |llj and |13| the authors consider a generalization of Lucas "Square Pyramid" 
problem, namely the determination of all pairs (a, b) £ Z^q that satisfy the equation 



b q 



l k + 2 k + 



+ a 



for some q > 2 and k > 1. Theorem 14.51 below is already proved in [2J where the authors determine 
all the solutions for all values of q and for 1 < k < 11. It is suggested in [T] that this could have 
been extended to larger values of k. Nevertheless, we present here this special case since our proof, 
which involves descent, includes the determination of the full set of rational points (as opposed to the 
subset of integral points) of a singular superelliptic curve of genus 7. We prove the case with q = 3 
because current tools only allow the use of cubic curves for the intermediate steps (although recent 
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developments, like the implementation of descent on Jq in MAGMA, suggest that Chabauty's method 
for superelliptic curves should be available in the near future, which would open the possibility to solve 
cases with q > 3). We consider the case with k = 9 since this is a value of k where the superelliptic 
curve involved in the computation is both singular and non-hyperelliptic. 

Theorem 4.5. The only pair (a,b) G Z 2 satisfying 

b 3 = ^i 9 
i=i 

IS (1,1). 

Proof. After replacing the right hand side of the equation by the closed formula for the sum of the first 
a ninth powers we get the equation 

b 3 = _L a 2 (a + 1)2(a 2 + a _ 1} ^4 + 2a 3 _ l fl 2 _ 3 fl+ 3\ (4 5) 

After the change of variables in the proof of Proposition 12. 3| we see that the solutions (a, b) correspond 
to rational points on the following genus 7, singular, superelliptic curve C in P 2 (l,4, 1) 

C : Y 3 = X 2 (X+5Z) 2 (X+WZ) 2 (X 2 +30XZ+W0Z 2 )(X 4 +30X 3 Z+4Q0X 2 Z 2 + 2A00XZ 3 +m0Z 4 ) 

(a, b) i-> (10, 100006, (a - 1)) G C(Q) (4.6) 
We have the descent map (jlq : C(Q) — > Aq/Q*Aq where the algebra Aq is isomorphic to the product 
QxQxQx K 1 xK 2 where K x = Q[t]/(t 2 + 30t + 100) and K 2 = Q[t]/(t 4 + 30t 3 + 460t 2 + 2400t + 4000). 
Denote the images of t in K\ and K 2 by 0\ and 2 respectively. We can compute representatives in 
for the images of the five known rational points. These are shown in Table 21 



i 


Pi G C(Q) 


m( p i) 


1 


(1,1,0) 


[(1,1,1,1,1)] 


2 


(0,0,1) 


[(l,5,10,ei,e 2 )] 


3 


(-5,0,1) 


(1, 75, 1, 1, 400(56*^ + 1346>! + 186O0 2 - 5840)) 


4 


(-10,0,1) 


(2, 1, 300, 25(30i + 10), 20(6^ + 200^ - 9O0 2 - 700)) 


5 


(-10,1000,3) 


(2, 1,4, 25(30i + 10), 20(130^ + 3300^ + 438O0 2 + 9600)) 



TABLE 4. Representatives in Aq of images of known points. 



By Proposition I3.12[ the genus G of the covers D a is 

G = q d ~ 2 ^ q ~ 1) - g )+l = 3 7 (^-3) +1 = 13123 

so the set S up (see Definition 13. 13j) contains the rational primes which are less than 4 x 13123 2 = 
688852516. Performing the local computations for all of these primes would be unfeasible, but we do 
not need to, since after checking the primes 2, 3 and 5 we exclude every element from Hq(S) apart 
from the five we already know. Thus 

SelW(C, Q) = hq ({Pi, P 2 , P 3 , P 4 , P 5 }) . 

These elements correspond to five covers 4>i : Di — )• C where each Di is a curve in P 8 x C, defined 
using a representative in /xq(Pj) (see (|3.2p in Section [3]). 

Now for each 1 < i < 5, we choose a combination of three factors from to form a cover 
Hi : Di — >• Ei where Ei is a genus one curve. By choosing factors whose product is defined over Q, we 
ensure that Ei and Kj are defined over Q. When choosing the three factors we also aim to obtain an 
elliptic curve Ei which has finitely many rational points. As it turns out, for all 1 < % < 5, the subset 
{x, x + 5, x + 10} C satisfies the criteria we need, giving the five curves in Table [5j 
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i 


(U, V, W)£Ei^ 


«i [i,Uh)her q ,(X,Y,Z)J = 


1 


V 3 = U{U + 5W)(U + WW) 




^ u fx)> Xu {x+5) U (x+lO) > Zu fx) J 




2 


50V 3 = U{U + 5W)(U + WW) 


\ Xu %) ' Xu Or+5)^(x+10) , z uf x) ^ 




3 


75U 3 = U(U + 5W)(U + WW) 




y Xu \x)^ Xu (x+h) u {x+W)-, Zu \ x ) j 




4 


mov 3 = u{u + 5W)(u + ww) 


\ Xu lx) ' Xu (x+5) U (x+W) ' ^ U (o;) > 




5 


8V 3 = U{U + 5W)(U + 10W) 


{ Xu fx)> Xu (x+5) u (x+l0), Zu \x) l 





Table 5. The cubic curves Ei and the maps «j : Dj — )■ E{. 



We note that there are isomorphisms E\ = E§ and E% = E^ over Q, but this is not relevant to the 
computation. We leave the coefficients of V 3 in E^ and E§ unchanged to remind the reader that they 
originate from multiplying the first three entries of the representatives of /Uq(Pj) shown in Table 21 
All five cubic curves are elliptic curves with Mordell-Weil rank equal to so we can determine the sets 



Ex( 

E 2 ( 
E 3 ( 
E4 

Ed 



(0,0,1), (-10,-10,3), (-20,10,3), 
10, 0,1), | (1 , 1 To)] , (-5,0, 1) 

-10,0,1), (-5,0,1)} 

]} 




-10,0,1), (-5,0, 1) 

(-10,0, 1) |, (-5,0,1)) 

(0,0,1), (-10, 0,1), (-5, 0,1), 
(-20, 5, 3), (2, 1,0), I (-10 ,-5 ,3) 



We then compute the pre-images of these sets under the maps Kj. The box indicates that a point has 
one Q-rational point in its Kj-fiber. The rest of the points have no rational pre-image. As expected, 



contains exactly one element for each 1 < i < 5. Using Corollary 

5 

= \J<t>i (D i (Q)) = {P 1 ,P 2 ,P 3 ,P i ,P 5 } 



we get 



i=l 

These points correspond to all solutions (a, b) G Q 2 satisfying Equation (|4.5p . In particular, using the 
inverse of the map (j4~6]) we get that Pi, P 2 , P3, P4 and P5 correspond to (1, 1), "00", (—1,0), (0,0) and 
(—2, 1) respectively. Therefore (1, 1) is the only solution where both a and b are positive integers. □ 
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